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Far-field  solutions  (source  mechanism  and  focal  depth)  have 
been  obtained  for  nineteen  events  in  the  Sinkiang,  China,  earthquake  swarm 
of  August  1974.  Both  spectral  and  ratio -of -events  fitting  methods  were  em- 
ployed. The  solutions  indicated  two  kinds  of  source  mechanisms  in  this 
swarm:  vertical  strike-slip  and  dip-slip  with  a dip  angle  of  forty  degrees. 

The  results  suggested  that  these  events  occurred  at  shallow  depths  with  the 
estimated  source  depth  varying  between  4 km  and  14  km. 

The  effective  Rayleigh  wave  attenuation  coefficients  along  sev- 
eral specific  travel  paths  have  been  determined  for  the  periods  of  10  seconds 
to  50  seconds.  The  two-station  method  was  applied  to  two  sets  of  observed 
data:  the  1974  Sinkiang  earthquake  swarm  and  earthquakes  used  in  the  previ- 
ous M -m.  study  (Turnbull  et  al.  , 1975).  This  study  yielded  results  for  the 
s b 

travel  paths  KON-OGD,  KON-ALQ,  KON-TLO,  and  TLO-ZLP  (where  KON, 
OGD,  TLO,  and  ZLP  are  VLPE  stations). 

The  original  and  a modified  version  of  the  Tryggvason  isotropic- 
source  method  for  the  determination  of  the  surface  wave  attenuation  has  been 
formulated  and  implemented.  They  were  employed  to  calculate  the  average 
Rayleigh  wave  attenuation  coefficients  along  the  travel  paths  between  the  east- 
ern Kazakh  and  three  array  stations  (ALPA,  LASA,  and  NORSAR),  and  those 
between  the  U.  S.  Nevada  test  site  (NTS)  and  the  three  array  stations.  All 
available  surface  wave  data  of  1975  EKZ  and  NTS  events  were  used.  The  re- 
sults indicated  that  the  attenuation  across  the  Eurasian  continent  was  much 
higher  than  that  across  the  North  American  continent  for  the  period  between 


20  and  50  seconds.  The  general  trend  of  the  attenuation  curve  obtained  here 
and  that  of  Tryggvason's  is  quite  similar,  although  the  former  shows  a slightly 
higher  attenuation. 
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SECTION  I 
INTRODUCTION 


During  this  report  period,  several  investigations  of  the  seismic 
source  from  far -field  surface  wave  data  have  been  continued.  For  the  source 
mechanism  study,  the  surface  waves  of  nineteen  events  in  the  Sinkiang,  China, 
earthquake  swarm  of  August  1974  are  analyzed,  and  the  results  are  discussed 
in  Section  II.  In  Section  III,  travel  path  effects  are  examined  where  several 
approaches  to  determine  the  surface  wave  energy  attenuation  are  presented. 
Using  the  surface  wave  data  of  earthquakes  (1974  Sinkiang  earthquake  swarm 
and  earthquakes  used  in  previous  Mg-m^  studies)  and  presumed  underground 
nuclear  explosions  in  1975  from  the  East  Kazakh  and  the  Nevada  test  site, 
Rayleigh  wave  energy  attenuation  coefficients  are  calculated  for  various 
travel  paths  of  interest. 
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SECTION  II 

ANALYSIS  OF  THE  SINKIANG,  CHINA,  EARTHQUAKE  SWARM 

OF  AUGUST  1974 

A.  INTRODUCTION 

Several  Eurasian  earthquakes  were  analyzed  and  discussed  in 
previous  reports  (Turnbull  et  al.  , 1973,  1974).  These  events  took  place  in 
Central  Eurasia  near  the  Chinese-Russian  border  and  were  quite  widely  sep- 
arated in  terms  of  their  locations  and  dates.  This  separation  gave  us  the 
chance  to  model  several  different  source  region  structures  and  to  study  the 
source  mechanisms  in  them. 

During  most  of  the  month  of  August  and  the  early  part  of  Sep- 
tember 1974,  a major  earthquake  swarm  of  over  300  recorded  events  occurred 
in  Sinkiang,  China.  The  geographic  location  of  the  swarm  is  shown  in  Figure 
II - 1 , with  the  travel  paths  encountered  to  the  VLPE  stations:  ALPA,  NORSAR, 
and  LASA.  Its  location  is  shown  again  by  the  boxed  area  in  Figure  II-2,  but 
this  time  in  relation  to  the  major  geologic  features.  A study  of  this  earth- 
quake swarm  offered  the  following  advantages: 

• The  opportunity  to  analyze  a large  number  of  events  from  a 
small  region,  investigating  the  repeatability  of  source  mechan- 
ism from  event  to  event. 

• The  possibility  of  using  a master  event  or  events  with  spectral 
ratios  to  eliminate  the  travel  path  effects,  which  usually  are 
complicated  and  not  well  understood;  and  the  station  instrument 
response  corrections,  which  sometimes  are  not  too  reliable 
(especially  those  of  VLPE  stations),  for  a better  estimation  of 
the  source  parameters. 
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LOCATION  OF  1974  SINK IANG  CHINA  EARTHQUAKE  SWARM 
RELATIVE  TO  THE  OBSERVATION  STATIONS 


LOCATION  OF  1974  SINKIANG  CHINA  EARTHQUAKE  SWARM 


• The  opportunity  to  determine  the  effective  surface  wave  atten- 
uation along  some  specific  travel  paths,  using  the  two -station 
method  with  more  than  one  event. 

In  this  section  are  presented  discussions  and  results  on  the 
source  parameter  estimation  for  the  chosen  events  in  this  swarm.  The  travel 
path  attenuation  using  the  available  data  from  the  swarm  will  be  given  in  the 
next  section  where  all  discussions  and  results  are  devoted  to  the  travel  path 
attenuation  with  other  additional  data  and  methods. 

B.  PROCEDURE  OF  ANALYSIS 

Nineteen  events  from  this  earthquake  swarm  have  been  chosen 
to  be  studied  here.  In  selecting  these  nineteen  events  out  of  about  sixty  pro- 
cessed events  for  this  swarm,  the  following  criteria  were  used: 

• The  event  is  a single  event  at  all  recording  stations;  that  is, 
the  mixed  events  due  to  the  overlapping  signal  time  gates  are 
excluded. 

• The  observed  data  are  available  at  three  or  more  stations  with 
good  data  quality.  Hopefully,  these  stations  are  distributed  in 
two  consecutive  quadrants. 

Data  quality  is  judged  by  the  appearances  of  the  time  traces,  the  dispersion 
curves,  and  the  uncorrected  amplitude  spectra.  The  typical  Rayleigh  wave 
time  traces,  dispersion  curves,  and  uncorrected  amplitude  spectra  are  given 
in  Appendix  A for  the  event  LX+SINK+SO 18. 

To  show  their  relative  locations,  Figure  II- 3 displays  the  de- 
tailed locations  of  these  events.  Their  locations,  dates,  origin  times,  body- 
wave  magnitudes,  and  PDE  depths  (if  available)  are  given  in  Table  II-l.  In- 
formation concerning  the  observation  stations  is  given  in  Table  II-2  for  each 
event. 
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FIGURE  II- 3 

DETAILED  LOCATIONS  OF  THE  EVENTS 
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T ABJ-.E  II- 1 

1974  SINKIANG  EARTHQUAKE  SWARM 


Origin 

Time 


Depth 


9 


Azimuth  at  source  ( + for  NE)  4 = Distance  to  source  in  degrees. 


OBSERVATION  STATIONS  FOR  1974  S1NKIANG  EARTHQUAKE  SWARM 

(PAGE  2 OF  2) 


Azimuth  at  source  (+  for  NE)  A = Distance  to  source  in  degrees. 


These  nineteen  events  are  discussed  in  two  groups,  as  indi- 
cated in  both  Tables  II- 1 and  II-2.  The  first  group  consists  of  five  events  for 
which  both  Rayleigh  and  Love  wave  data  are  available.  The  number  of  obser- 
vation stations  is  betweeen  six  and  eight  for  this  group.  The  remaining  four- 
teen events  constitute  the  second  group.  It  is  noticed  that  observation  stations 
are  distributed  in  the  NE  and  the  NW  quadrants  for  all  nineteen  events.  The 
procedure  used  in  the  analysis  of  the  events  is  as  follows: 

• Choice  of  source  region  earth  model  - The  layered  earth  model 
is  determined  by  comparing  the  source  region  earth  structure 
profiles  with  the  known  earth  models  (as  discussed  by  Turnbull 
et  al.  , 1974).  The  Harkrider's  medium  response  solutions  are 
calculated  for  the  selected  earth  model  in  order  to  produce  the 
theoretical  amplitude  spectra.  For  this  swarm,  it  has  been 
found  that  the  normal  Gutenberg-Bullen  earth  model  is  most 
appropriate. 

• Data  handling  - Applying  a series  of  narrowband  filters  (as 
discussed  by  Turnbull  et  al.  , 1974),  each  component  is  ana- 
lyzed for  evidence  of  multipathing.  At  each  period,  group 
velocities  are  determined  for  each  multipath,  and  then  com- 
pared to  the  standard  group  velocity  curves  for  the  type  of 
travel  path  encountered  (see  example  in  Appendix  A).  Estimat- 
ing the  correct  multipath  from  this  comparison,  the  amplitude 
spectra  is  obtained  and  then  corrected  for  attenuation,  geomet- 
ric spreading,  and  instrument  response. 

• Source  parameter  determination  - Spectral  fitting  procedures 
(Turnbull  et  al.  , 1973)  are  then  applied  to  the  corrected  spec- 
tra, using  a theoretical  source  model  placed  in  the  selected 
layer  earth  model.  Source  parameters  are  estimated  by  first 
using  simple  spectral  fitting.  A master  (or  reference)  event 
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(or  events)  is  picked  based  on  the  information  obtained  from  the 
spectral  fitting.  Then  the  ratio-of-events  fitting  is  applied  to 
the  rest  of  the  events  in  the  same  group  using  the  selected  mas- 
ter event  (or  events)  as  reference. 

C.  FAR-FIELD  SOLUTION 
1.  Group-I 

For  this  group,  observed  surface  wave  data  are  available  for 
both  Rayleigh  and  Love  waves  at  most  of  the  stations.  The  observation  sta- 
tions include  all  three  array  stations:  ALPA,  LASA,  and  NORSAR.  In  gen- 
eral, data  quality  is  better  at  the  array  station  than  at  the  VLPE.  Here,  we 
first  apply  the  spectral  fitting  method  using  all  available  Rayleigh  and  Love 
wave  data.  For  the  five  events  in  this  group,  the  estimations  of  the  source 
parameters  based  on  the  minimum-residual  criterion  are  listed  in  Table  II-3, 
and  those  based  on  the  distribution-of-minimum- residual  criterion  are  listed 
in  Table  II-4.  The  distributions  of  minimum  residuals  of  each  parameter  of 
each  event  are  shown  in  Figures  II-4  through  U-8  with  their  spectral  fits. 

Referring  to  those  figures  of  the  spectral  fits,  it  is  noticed  that 
the  spectral  fits  for  the  Love  wave  in  general  are  poor  in  both  spectral  level 
and  in  shape,  except  for  the  events  LX+SINK+S050  and  S052.  This  is  probably 
caused  by  the  poor  quality  of  the  observed  Love  wave  spectra.  To  determine 
whether  the  estimation  of  the  source  parameters  will  be  significantly  affected 
by  these  poor  Love  wave  data,  we  next  apply  the  spectral  fitting  method  using 
Rayleigh  wave  data  only.  The  results  from  this  fitting  are  given  in  Tables 
II - 5 and  II - 6 and  in  Figures  II - 9 through  II - 1 3 . 

Comparing  the  results  obtained  from  two  spectral  fittings  (see 
Tables  II - 3 and  II - 5) , using  both  Rayleigh  and  Love  wave  and  using  Rayleigh 
wave  only,  no  significant  differences  in  the  source  parameter  estimations 
have  been  observed,  except  the  dip  and  slip  angles  for  events  LX+SINK+S001 
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ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED 
BY  AMPLITUDE  SPECTRAL  FITTING  BASED 
ON  MINIMUM-RESIDUAL  CRITERION: 

ALL  AVAILABLE  LR  AND  LQ  DATA 
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Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

x° 

Strike 

N^°E 

Moment 

..25, 

10  dyne-cm 

LX+SINK+S001 

6 

40 

50 

+90 

75 

0.  123  x 102 

LX+SINK+S018 

6 

60 

60 

80 

0.  161  x 101 

LX  +SINK+S024 

6 

40 

-60 

130 

0.  193 

LX+SINK+S050 

4 

40 

50 

±90 

35 

0.  797 

LX+SINK+S052 

4 

40 

50 

+90 

60 

0.  827  x 10"1 
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SOURCE  PARAMETER  DISTRIBUTIONS: 
LX+SINK+S001  - BOTH  LR  AND  LQ 
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SPECTRAL  FIT:  LX+SINK+SOOl  - BOTH  LR  AND  LQ 
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SPECTRAL  FIT:  LX+SINK+S018  - BOTH  LR  AND  LQ 
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TABLE  II - 5 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED 
BY  AMPLITUDE  SPECTRAL  FITTING  BASED 
ON  MINIMUM-RESIDUAL  CRITERION: 

LR  ONLY 


i 


Event  1.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

x° 

Strike 

N0°E 

Mome  nt 
10  dyne-cm 

LX+SINK+S001 

10 

90 

0 

m 

0.  188  x 10Z 

LX+S1NK+S018 

6 

50 

-30 

— 

120 

0.  198  x 101 

LX  +SINK+S024 

10 

90 

0 

40 
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0.  597 

LX  +SINK+S050 

4 

40 

50 

±90 

55 

0.  707 

LX+SINK+S052 

4 

40 

50 
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55 
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ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE  SPECTRAL 
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and  S024.  In  fact,  for  events  LX+  SINK+  SO 50  and  SO 52,  the  Love  wave  fits  of 
which  are  compatible  with  their  Rayleigh  wave  fits,  two  spectral  fittings  have 
yielded  practically  the  same  source  parameter  estimation.  For  events  LX  + 
SINK+S001  and  S024,  the  Love  wave  fits  of  which  are  poor,  the  spectral  fit- 
ting using  Rayleigh  wave  only  has  changed  the  source  mechanism  solution 
from  dip-slip  to  vertical  strike-slip.  This  difference  can  be  reasonably  ex- 
plained based  on  the  observed  fact  that  the  theoretical  Rayleigh  wave  ampli- 
tude spectra  of  the  pure  dip-slip  (slip  angle  of  90°)  and  the  vertical  strike-slip 
fault  are  very  similar  in  shape,  and  the  addition  of  the  poor  Love  wave  data  in 
the  spectral  fitting  does  not  help  in  resolving  this  similarity.  Although  based 
on  the  minimum-residual  criterion,  two  spectral  fittings  have  resulted  in 
quite  different  source  mechanisms  for  two  events  in  this  group.  Based  on  the 
distribution-of-minimum-residual  criterion,  these  fittings  give  more  or  less 
the  same  probable  ranges  in  the  source  parameter  estimation. 

Among  these  five  events,  events  LX+SINK+S018  and  S052  are 
chosen  to  be  used  as  the  reference  event  in  the  ratio-of-events  fitting,  based 
on  the  observation  that  the  spatial  fitting  solutions  for  these  two  events  seem 
to  be  more  stable  and  yield  a better  spectral  fit.  The  ratio-of-events  fitting 
method  is  applied  to  the  event  LX+SINK+S001  using  LX+SINK+S018  as  a refer- 
ence and  to  events  LX+SINK+S024  and  S050  using  LX+SINK+S052  as  a refer- 
ence. Referring  to  Table  II-2,  these  pairs  of  events  are  picked  such  that  two 
events  in  the  same  pair  have  the  maximum  number  of  common  observation 
stations.  Estimations  of  the  source  parameters  obtained  by  the  ratio-of- 
events  fitting  are  presented  in  Tables  U-7  through  II-9,  with  the  distributions 
of  minimum-residuals  shown  in  Figures  11-14  through  11-16.  For  the  event 
LX+ SINK+ SO  50,  the  ratio-of-events  fitting  has  yielded  an  estimation  of  the 
same  source  parameter  as  the  spectral  fitting.  For  the  event  LX+ SINK+ S00 1 , 
the  source  parameter  estimation  by  the  ratio-of-events  fitting  favors  the 
strike-slip  fault,  while  for  the  event  LX+  SINK+  S024,  the  estimation  favors 
the  dip-slip  fault. 


■Mi 
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TABLE  11-7 


ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  R AT lO-OF- EVENTS  FITTING 

Reference  Event  : LX+SINK+S018 

Solutions  for  Event:  LX+S1NK+S001  (SK+01/R-18) 


A.  Solution  by  Minimum -Residual  Criterion 


Event  I.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

A° 

Strike 

N^°E 

Moment 
1025  dyne -cm 

LX  +SINK+S001 

8 

80 

-30 

30 

0.  148  x 102 

B.  Solution  by  Distribution-Of-Minimum-Residual  Criterion 


Event  I.  D.  : LX+SINK+S001 

Source  Parameters 

Probable  Range 

% Confidence 

h 

6-10 

84 

6 

70-90 

79 

K 

-30, 0, +30 

93 

+ 

25-40 

23 

115-130 

38 

TABLE  II - 8 


ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  RATIO-OF -EVENTS  FITTING 

Reference  Event  : LX+SINK+S052 

Solutions  for  Event:  LX+SINK+S024  (SK+24/R-52) 


A.  Solution  by  Minimum-Residual  Criterion 


Event  I.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

A? 

Strike 

N^°E 

Moment 
1025  dyne -cm 

JLX+SINK+S024 

4 

40 

50 

+90 

65 

0.  179 

B.  Solution  by  Distr ibution-Of-Minimum-Residual  Criterion 


Event  I.  D.  : LX+SINK+S024 

Source  Parameters 

Probable  Range 

% Confidence 

h 

0.  5-4 

90 

6 

30-60 

60 

X 

- 

- 

4> 

- 

- 
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TABLE  II - 9 


ESTIMATIONS  OF  SOURCE  PARAMETERS 
OBTAINED  BY  R ATIO-OF  -EVENTS  FITTING 

Reference  Event  : LX+SINK+S052 

Solutions  for  Event:  LX+S1NK+S050  (SK  + 50/R-52) 


A.  Solution  by  Minimum -Residual  Criterion 


Optimal  Solution 

Event  I.  D.  Depth  Dip  Angle 

h (km)  <5  ° 


LX+SINK+S050 


Slip  Angle 

A° 

Strike 

N0°E 

Moment 
1 02 **  dyne  -cm 

±90 

65 

0.  828 

Solution  by  Distr  lbution-Of  - Minimum -Res  idual  Criterion 


Event  1.  D.  : LX+S1NK+S050 


Source  Parameters  Probable  Range 


0.  5-6 


% Confidence 


84 


82 


8 


STRIKE 


FIGURE  11-14 


SOURCE  PARAMETER  DISTRIBUTIONS  BY  R ATIO-OF- EVENTS: 

LX  +SINK+S001 
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The  major  results  from  the  analysis  of  this  group  are  as 

follows: 

• The  fitting  using  Rayleigh  wave  only  and  that  using  both  Ray- 
leigh and  Love  waves  should  give  more  or  less  the  same  source 
parameter  estimation.  When  both  Rayleigh  and  Love  wave  data 
are  good  and  their  data  quality  are  compatible,  we  can  expect 
to  have  more  resolution  in  the  strike  direction  by  using  both 
Rayleigh  and  Love  wave.  This  is  thought  to  be  expected,  since 
the  shape  of  the  theoretical  amplitude  spectra  for  the  Rayleigh 
wave  has  more  variations  with  respect  to  the  source  parameters 
than  that  for  the  Love  wave  (Turnbull  et  al.  , 1973),  and  the 
theoretical  Love  wave  radiation  patterns  of  various  periods  of 
interest  (10  to  50  seconds)  and  have  sharper  nodes  for  the  bet- 
ter diagnosis  of  the  strike  direction.  However,  when  the  Love 
wave  data  are  poor  as  compared  to  the  Rayleigh,  the  addition  of 
the  Love  wave  data  in  the  spectral  fitting  sometimes  will  change 
the  estimation  of  the  dip  and  slip  angle  very  noticeably. 

• The  estimated  source  parameters  for  this  group  indicate  that 

five  events  all  occurred  at  shallow  depth  (4  km  - 10  km)  with 
two  possible  kinds  of  source  mechanisms:  dip-slip  and  vertical 

strike-slip  faults. 

2.  Group-II 

There  are  fourteen  events  in  this  group.  However,  for  this 
group,  observed  surface  wave  data  are  available  only  for  the  Rayleigh  wave, 
and  the  number  of  observation  stations  is  between  three  to  five  with  all  three 
array  stations  available  for  only  six  events.  The  data  quality  for  this  group, 
nevertheless,  is  as  good  as  that  of  the  Rayleigh  wave  for  the  first  group. 

• 

Similar  to  the  first  group,  we  first  apply  the  spectral  fitting 
method  to  each  event.  The  estimations  of  the  source  parameters  based  on 
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the  minimum-residual  criterion  are  listed  in  Table  II- 10,  and  those  based  on 
the  distribution -of -minimum -residual  criterion  are  listed  in  Table  U-ll.  The 
distributions  of  minimum  residuals  of  each  parameter  of  each  event  are  given 
in  Appendix  B.  These  results  again  indicate  that  for  this  group  there  are  two 
kinds  of  source  mechanisms  --dip-slip  and  vertical  strike-slip  fault  --  and 
the  focal  depths  of  these  events  are  all  shallow,  ranging  between  4 km  and  14 
km. 

We  also  apply  the  ratio -of-events  fitting  method  to  these  events 
using  the  event  LX+ SINK+ SO  10  as  a reference.  The  event  LX+ SINK+ S010  is 
chosen  to  be  the  reference  because:  (1)  this  event  has  better  spectral  fits, 

especially  for  three  array  stations;  (2)  the  estimation  of  the  focal  depth  from 
the  spectral  fitting  agrees  quite  well  with  the  PDE  depth;  and  (3)  all  other 
events  in  this  group  have  the  largest  number  of  common  observation  stations 
with  this  event.  The  results  obtained  from  the  ratio -of-events  fitting  method 
are  presented  in  Tables  11-12  and  11-13. 

Referring  to  Tables  II - 10  and  11-12,  which  are  based  on  the 
minimum-residual  criterion,  there  is  good  agreement  between  the  focal  depth 
estimation  from  the  spectral  fitting  and  that  from  the  ratio-of-events  fitting. 
Although  for  several  events  two  fittings  have  yielded  quite  different  estima- 
tions of  the  dip  and  the  slip  angle  based  on  the  minimum-residual  criterion, 
the  probable  ranges  of  the  dip  and  the  slip  angle  estimations  (see  Tables  11-11 
and  11-13)  obtained  from  the  fitting  are  very  similar.  This  is  the  same  situ- 
ation that  was  observed  and  discussed  in  the  first  group. 

In  Figure  11-17,  we  have  plotted  the  estimated  seismic  moment 

(m  ) versus  bodywave  magnitude  (m,  ) for  each  event,  in  both  the  first  and  the 

2 b 

second  group,  in  relation  to  the  U -model  determined  by  Tsai  (1972).  We 
see  that  there  is  close  agreement. 
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TABLE  II- 10 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED 
BY  AMPLITUDE  SPECTRAL  FITTING  BASED 
ON  MINIMUM-RESIDUAL  CRITERION 
(PACE  1 OF  2) 


Event  1.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

x° 

Strike 

N0°E 

Moment 

lrt25_, 

10  dyne-cm 

LX  +S1NK+S004 

10 

90 

0 

30 

120 

0.  391 

LX  +S1NK+S005 

4 

40 

-60 

100 

0.  1 17 

LX  +SINK+S006 

12 

90 

0 

30 

120 

0.  496 

LX  -&S1NK+S008 

6 

40 

-60 

95 

0.  224 

LX+SINK+S010 

10 

90 

0 

30 

120 

0.  195 

LX  +SINK+S01  1 

10 

90 

0 

25 

115 

0.  265  x 10'1 

LX+SINK+S017 

10 

90 

0 

30 

120 

0.  171  x 101 
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TABLE  II- 10 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED 
BY  AMPLITUDE  SPECTRAL  FITTING  BASED 
ON  MINIMUM-RESIDUAL  CRITERION 
(PAGE  2 OF  2) 


Event  I.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

6° 

Slip  Angle 

k° 

Strike 

N^°E 

Moment 
, 25 

10  dyne-cm 

LX  +SINK+S02  3 

8 

80 

0 

125 

0. 967  x 10'* 

LX+S1NK+S025 

4 

40 

50 

+90 

70 

0.  254  x 10'1 

LX+SINK+S031 

8 

90 

0 

30 

120 

0. 571  x 10'1 

LX+SINK+S036 

4 

40 

50 

+90 

65 

0.475  x 10_1 

LX+SINK+S04  3 

10 

80 

0 

120 

0.  423  x 10'1 

LX+SINK+S056 

14 

90 

0 

40 

130 

0.  387 

LX+SINK+S059 

14 

90 

0 

E9 

0.  147 

ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE  SPECTRAL 
FITTING  BASED  ON  DISTRIBUTION-OF -MINIMUM-RESIDUAL  CRITERION 

(PAGE  1 OF  2) 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  AMPLITUDE  SPECTRAL 
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TABLE  11-12 


ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED 
BY  R ATIO-OF -EVENT  FITTING  BASED  ON 
MINIMUM-RESIDUAL  CRITERION 
(PAGE  1 OF  2) 


Reference:  LX+SINK+S010  h = 1 Okm,  6 = 90,  A.=  0,  <f>=  30, M=  0.  1 95 


Event  I.  D. 

Optimal  Solution 

Depth 
h (km) 

Dip  Angle 

<5° 

Slip  Angle 

X° 

Strike 

N^°E 

Moment 

25J 

10  dyne-cm 

LX  +SINK+S004 

10 

90 

0 

H 

0.  426 

LX  +SINK+S005 

4 

m 

±90 

105 

0.  893  x 10'1 

LX+SINK+S006 

10 

90 

0 

25 

115 

0.  675 

LX  +S1NK+S008 

8 

80 

0 

85 

0.  262 

LX+SINK+S01  1 

6 

40 

50 

±90 

60 

0.  133  x 10_1 

LX+SINK+S017 

10 

90 

0 

ia 

0.  194  x 101 
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Refe  rence : 

LX+SINK 

Event  I.  D. 

Depth 
h (km) 

LX  +S1NK+S02  3 

6 

LX  +SINK+S02  5 

6 

LX  +SINK+S03 1 

6 

LX  +SINK+S036 

8 

LX  +S1NK+S043 

6 

LX  +SINK+S059 

6 

i 

i 


jL 


Optimal  Solution 


Dip  Angle 

6° 

Slip  Angle 

A° 

Strike 

N^°E 

Moment 

25 

10  dyne-cm 

80 

0 

35 

0.  676  x 10'1 

90 

0 

H 

0.  451  x 10'1 

80 

0 

30 

0.  318  x 10‘* 

90 

0 

25 

0.  177 

H 

±90 

60 

0.  258  x 10'1 

90 

0 

25 

115 

0.  454  x 10'1 
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ESTIMATIONS  OF  SOURCE  PARAMETERS  OBTAINED  BY  R ATIO-OF- EVENT 
FITTING  BASED  ON  DISTRIB  UTION-OF  - MINIMUM- RESIDU  AL  CRITER  ION 

(PAGE  1 OF  2) 
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D.  CONCLUSIONS 


Analysis  of  the  far-field  surface  wave  spectra  of  the  nineteen 
events  in  the  1974  Sinkiang,  China,  earthquake  swarm  has  yielded  consistently 
shallow  focal  depth  estimations  ranging  from  4 km  to  1 4 km  for  these  events. 
This  agrees  quite  well  with  the  PDE  depth  which  is  available  for  five  events 
(see  Table  II - 1 ) . Also,  from  the  distributions  of  the  minimum-residuals  for 
the  focal  depth  estimations  obtained,  it  is  believed  that  the  shallow-focal 
depth  estimations  should  be  quite  reliable.  For  the  estimations  of  the  dip  and 
the  slip  angles,  the  results  suggest  that  there  were  two  possible  kinds  of 
source  mechanisms  involved  in  this  earthquake  swarm:  The  dip-slip  fault 

probably  was  responsible  for  six  events  and  the  other  thirteen  events  had 
vertical  strike-slip  faults.  It  is  conjectured,  however,  that  the  vertical 
strike-slip  fault  is  the  more  probable  source  mechanism  for  events  in  this 
earthquake  swarm  based  on  the  following  observed  facts: 


• The  source  parameter  estimations  have  resulted  in  the  vertical 
strike-slip  fault  for  the  majority  of  the  events  (thirteen  out  of 
nineteen)  being  analyzed  for  this  earthquake  swarm.  More- 
over, for  these  events  with  the  dip-slip  fault,  the  distributions 
of  the  minimum-residuals  indicate  that  the  strike-slip  fault  is 
also  possible. 

• The  theoretical  amplitude  spectra  of  the  vertical  strike-slip 
fault  and  those  of  the  normal  dip-slip  fault  (i.  e.  , slip  angle  of 
+ 90°)  are  very  similar  in  shape. 

• The  previously-analyzed  earthquakes  in  this  region  seem  to 
favor  the  strike-slip  fault  (Tsai  and  Shen,  1972;  Turnbull  et 
al.  , 1973,  1974). 


The  estimation  of  the  strike  direction  for  the  events  with  vertical  strike-slip 
faults  is  about  N40°~50°W  (or  N50°~40°E).  This  is  quite  consistent  from 
event  to  event. 
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SECTION  III 

TRAVEI,  PATH  ATTENUATION 


A.  INTRODUCTION 

There  are  several  factors  which  can  affect  the  spectral  level 
and  shape  of  the  observed  far-field  surface  wave  amplitude  spectra  of  a given 
seismic  event.  These  factors  can  be  conveniently  lumped  into  three  groups: 
the  source  mechanism,  the  travel  path  effect,  and  the  instrument  response  of 
the  observation  station.  Therefore,  in  order  to  estimate  the  seismic  source 
parameters  using  the  observed  far-field  surface  wave  spectra,  some  appropri- 
ate corrections  to  the  observed  spectra  must  be  made.  The  corrections  for 
the  instrument  response  and  the  geometric  spreading,  which  are  considered 
part  of  the  travel  path  effect  here,  are  trivial;  while  the  correction  for  the 
energy  attenuation,  which  is  the  other  part  of  the  travel  path  effect,  generally 
is  not  well  known  due  to  the  limited  knowledge  of  the  surface  wave  energy  dis- 
sipation. 

The  travel  path  attenuation  correction  is  important  in  the  esti- 
mation of  the  seismic  source  parameters  by  the  spectral  fitting  method  using 
observed  far-field  surface  wave  data.  The  level  of  the  attenuation  curve  as 
the  function  of  frequency  can  affect  the  estimation  of  the  seismic  moment,  and 
the  general  shape  of  the  curve  will  affect  the  estimation  of  the  focal  depth  and 
source  mechanism.  So  far,  when  no  reliable  knowledge  of  energy  attenuation 
along  the  travel  paths  of  interest  is  available,  the  Tryggvason's  energy  atten- 
uation curve  for  the  normal  continental  path  (Tryggvason,  1965)  has  been  used. 
The  general  shape  of  the  Tryggvaron  curve  is  believed  to  be  reasonable  while 
the  level  of  the  curve  seems  too  low  (Tryggvason,  1965;  Anderson,  1964).  It 
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is  the  purpose  of  this  section  to  obtain  the  surface  wave  attenuation  curves 
for  several  travel  paths  of  interest. 

Several  methods  to  determine  the  surface  wave  attenuation  coef- 
ficient will  be  derived  and  discussed  in  Subsection  B.  In  Subsection  C,  these 
methods  are  applied  to  several  groups  of  observed  surface  wave  data,  and 
Rayleigh  wave  energy  attenuation  curves  are  obtained  for  several  travel  paths 
of  interest. 


B.  METHODS 

The  term  'energy  dissipation'  or  'energy  attenuation'  used 
throughout  this  section  is  chosen  to  indicate  all  energy  lost  along  the  travel 
path  including  absorption,  reflected  energy,  and  scattered  energy. 

In  the  following  derivation,  it  is  assumed  that  there  is  no 
energy  transferred  from  one  frequency  to  the  other  along  the  travel  path  and 
that  the  observed  amplitude  spectrum  has  been  corrected  for  the  station  in- 
strument response.  Two  methods  will  be  described  below:  the  two-station 

method  which  can  be  employed  by  the  earthquake  data,  and  the  isotropic - 
source  method  which  is  suitable  for  the  explosion  data. 

1.  Two -Station  Method 


The  observed  surface  wave  amplitude  spectra  of  a given  seis- 
mic event  recorded  at  two  stations  can  be  related  as  follows: 


Aj  (f ) 

vd 


R(oi,o2,f) 


/ \ 1 /2 

/sinA^V  e-kA1(f)D1+kA2(f)D2 
l sinAj  I 


(HI-  1 ) 


where 

A^(f)  is  the  observed  amplitude  spectrum  for  frequency  f at  station  i, 
k^.(f)  is  the  amplitude  attenuation  coefficient  for  frequency  f along 
the  travel  path  i (great  circle  path  oetween  the  seismic  event 
and  station  i ), 
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D. 

1 


O. 

1 


R^°r  °2’ f) 


is  the  epicenter  distance  in  km  to  station  i , 
is  the  epicenter  distance  in  degrees  to  station  i , 
is  the  source  azimuthal  angle  to  station  i , and 
is  the  ratio  of  the  source  excitation  in  to  that  in 

for  frequency  f . 


In  equation  (HI-  1 ) the  term  R(0j,0.,,f)  accounts  for  the  non- 
circular radiation  pattern  of  the  source.  For  an  isotropic  source,  R(0j,0,  , f) 
is  one.  For  an  anisotropic  source,  such  as  an  earthquake,  in  general, 

R(0  , O , f)  is  different  from  one.  However,  R(0  , O , f)  will  be  one  if  O 
is  equal  to  , regardless  of  source  type.  That  is,  when  two  observation 
stations  are  on  the  same  great  circle  path  between  the  seismic  event  and  the 
station  and  have  the  same  azimuthal  angle,  equation  ( III-  1 ) can  be  reduced 
as  follows: 


*,<<>  -kA|,HDl-D2> 

A2(f)  " 


(111-2) 


along  the  travel  path  between  stations  1 and  2.  Therefore,  equation  ( 111-2 ) 
can  be  applied  to  the  appropriate  earthquake  data  to  calculate  the  travel  path 
attenuation  as  follows: 


k£(f)  = 2kA(f)  = 


(III-  3) 


where  k^,(f) 
E 


is  the  energy  attenuation  coefficient  for  frequency  f . 
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2. 


Isotropic  Source  Method 


Tryggvason’s  approach  assumes  that  radiation  is  isotropic  and 
that  attenuation  is  the  same  along  all  event-station  travel  paths.  Then  the 
energy  attenuation  coefficient  can  be  estimated  from  the  following  equation  by 
using  the  far-field  observed  surface  wave  spectra  at  several  observation  sta- 
tions for  one  explosion  (Tryggvason,  1965): 


2 f n A. (f)  + in  sinA.  = C^i)  - D.k^f)  i = 1 , NSITE  (III-4) 

where  A.,  A.,  D.,  and  k£  have  the  same  meaning  as  before.  C ^ is  a con- 
stant depending  on  frequency  and  the  size  of  explosion.  NSITE  is  the  number 
of  observation  stations.  In  equation  (1II-4),  k and  C,  are  two  unknown 
parameters  to  be  estimated.  Hence,  when  NSITE  is  greater  than  two,  k 
and  C2  can  be  estimated  by  the  least-squares  method,  i.  e.  , minimizing  the 
following  quantity: 

NSITE 

€ = [(C2  " °ikE)  ' (2^nA.  + /nsinA.)J2  . (III-5) 


The  estimated  k will  be  as  follows: 
E 


NSITE  Sjp/nA.  + /n  sin  A.)dJ-(£  D.)[S  (ZtnA.  + in  sin  A.)] 

D,)(Ld.)  - NSITE  £ (D  D ) 
iiii  i i i' 


(III  - 6 ) 


whe  re 


E 


stands  for 


NSITE 

E • 


i = l 


To  facilitate  the  later  discussion,  this  original  Tryggvason's  approach  will  be 
named  as  ISM-A.  It  can  be  easily  shown  that  equation  (1II-6)  will  be  reduced 
to  equation  (111-3)  of  the  two-station  method  when  NSITE  is  equal  to  two,  as 
expected. 
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The  ISM-A  utilizes  only  the  observed  surface  wave  spectra 

from  one  explosion.  When  there  are  several  explosions  taking  place  in  the 

same  area,  the  best  usage  of  these  available  data  by  the  ISM-A  will  be  that 

when  k is  estimated  from  equation  (III-6)  by  using  the  observed  spectra  of 

one  explosion  j , and  the  final  estimate  of  k is  obtained  by  some  curve- 

E 7 

fitting  technique  or  just  by  taking  the  following  simple  average: 


NEXP 

kE  ' £ VNEXP  • 


(III-7) 


where  NEXP  is  the  number  of  available  explosions.  When  there  is  a suffic- 
iently large  number  of  the  observation  stations,  the  variation  between  k 's 

E 

and  kEj'S  Can  be  exPected  to  be  small.  However,  when  the  number  of  obser- 
vation stations  is  too  small,  say  three,  large  variations  among  k 's  can  be 

Ej 

expected.  For  this  case,  to  take  the  advantage  of  several  available  explosives 
and  to  best  utilize  the  observed  data,  we  can  modify  the  ISM-A  by  simultane- 
ously taking  all  explosions  into  account  in  the  estimation  of  k as  follows. 

E 

Modification  of  equation  ( III-4 ) by  considering  all  available  ex- 
plosions at  the  same  time  yields  the  following  result: 


2fn  A..  + fnsinA..  - C_  . - D k 
ij  ij  2j  ij  E 


i = 1 , NSITE  j = 1 , NEXP  ( III  - 8 ) 

where  now: 


A. 


is  the  observed  amplitude  spectrum  for  frequency  f at  sta- 
tion i due  to  the  explosion  j , 

is  the  epicenter  distance  in  km  and  degrees,  respectively, 
between  the  station  i and  the  explosion  j , 
is  a constant  for  frequency  f and  explosion  j. 
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ing  quantity: 


Proceed  as  in  the  ISM- A,  but  this  time  minimize  the  follow- 


€ 


NEXP  NSITE 

Z Z [(C2j'DijkE)'{2/nAij  + i'nsiriAij)J  * (III'9) 


The  estimated  k will  be  as  follows: 


NSITE  ZECzf  nA. . +/nsinA  ,)(D.  )]  - SnC(2  fnA  + fnsinA  )1(EZd  ) 

J * l iJ U ij  J J Li ij ij  J i ij; 

£["(Sd..)(£d..)  1 - NSITE  ee(DD) 
j L i ij  i ij  J j i ij  ij' 


where 


stand  for 


NEXP 

and 

j = l 


NSITE 


£ 


i=l 


respec  tively. 
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We  will  call  this  modified  approach  an  ISM-B.  With  ISM-B, 
as  long  as  NSITE  ^ 2 and  NEXP  >2,  the  energy  attenuation  coefficient  k 

£ 

can  be  estimated;  while  with  ISM-A,  NSITE  must  be  greater  than  two  in  order 

to  estimate  k^..  When  NSITE  = 2,  the  ISM-A  can  offer  only  the  deterministic 

value  of  k„  , as  does  the  two-station  method.  With  a small  number  of  avail- 
E 

able  observations  stations,  therefore,  the  ISM-B  probably  will  yield  a better 
estimate  of  k , as  long  as  NEXP>  2. 


C.  RESULTS 

To  calculate  the  surface  wave  energy  attenuation  for  the  travel 
paths  of  interest,  the  methods  discussed  in  the  previous  subsection  are  ap- 
plied to  two  classes  of  observed  data:  Rayleigh  waves  of  the  earthquake  and 
of  the  explosion.  Here,  before  the  application  of  these  methods,  the  observed 
Rayleigh  wave  data  are  preprocessed  by  narrowband  filtering  to  eliminate 
any  possible  muitipathing,  as  done  in  Section  11.  The  observed  Rayleigh  wave 
spectra  are  obtained  and  then  corrected  for  the  station  instrument  response. 
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1. 


Earthquake  Data 


Two  sets  of  earthquake  data  --  the  1974  Sinkiang,  China, 

earthquake  swarm  and  earthquakes  used  in  the  previous  M -m,  study  (Turn- 

s b 

bull  et  al.  , 197  5)  --  were  examined  to  determine  which  events  could  be  em- 
ployed by  the  two-station  method.  As  previously  mentioned  (see  Subsection 
B-l),  with  earthquake  data,  the  two -station  method  can  be  applied  to  the  events 
which  have  two  observation  stations  lying  on  the  same  great  circle  path  with 
the  same  source  azimuth. 

For  the  1974  Sinkiang  swarm,  there  were  five  events  which 
approximately  satisfy  the  above  requirement.  The  pertinent  information  about 
the  events  and  stations  for  these  five  events  is  given  in  Table  III  - 1 with  the 
event-station  travel  paths  plotted  in  Figure  III  - 1 . The  two  stations  involved 
are  VLPE  stations  at  Toledo,  Spain  (TLO),and  La  Paz,  Bolivia  (ZLP).  It  is 
noticed  that, although  these  two  stations  are  not  quite  on  the  same  event-station 
great  circle  path,  their  event-station  travel  paths  are  fairly  close:  two  source 

azimuths  differing  by  3.  5 degrees.  The  calculated  energy  attenuation  coeffic- 
ients along  this  travel  path  (TLO-ZLP)  are  plotted  in  Figure  III-2,  with  the 
variation  among  events  given  as  a vertical  bar  and  the  average  values  of  k 

E 

for  various  frequencies  connected  by  a solid  line. 

For  earthquakes  used  in  the  previous  Mg-m^  study  (59  events), 
three  events  were  found  suitable  for  the  two-station  method.  The  pertinent 
information  about  the  events  and  the  stations  is  listed  in  Table  III-2  with  the 
event-station  travel  paths  plotted  in  Figure  III- 3.  Three  travel  paths  for 
which  the  two -station  method  can  be  used  to  calculate  the  travel  path  attenu- 
ation are  as  follows: 

• KON  (VLPE  station  at  Kongsberg,  Norway)  to  TLO  obtainable 
from  the  event  LX+ KURSP+ 425, 

• KON  to  OGD  (VLPE  station  at  Ogdensburg,  New  Jersey)  ob- 
tainable from  the  event  LX  + CAUCA  + 1 24,  and 
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TABLE  III- 
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EVENT-STATION  TRAVEL  PATH:  TLO-ZLP 
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TABLE  III-2 

INFORMATION  OF  EVENTS  AND  STATIONS 
(a)  Path  KON-TLO 


Event:  LX+KURSP+425  ; 04/05/73  ; 43.  6N  , 147.  7E  ; mb=  5.4 


Station 

Latitude 

°N 

Longitude 

°E 

* 

AZ 

** 

4° 

A km 

KON 

59.  7 

9.  59 

-20.  9 

71.  3 

7937.  9 

TLO 

39.  7 

1 

4^ 

O 

ro 

-21.3 

92.  9 

10336. 2 

(b)  Path  KON-OGD 


Event:  LX-CAUCA-124  ; 02/03/72  ; 40.  ON  , 48.  4E  ; mb=  5.  1 

Station 

Latitude 

°N 

Longitude 

°E 

* 

AZ 

A° 

A km 

KON 

59.  7 

9.  59 

-37.  9 

31.  1 

3464. 4 

OGD 

41.  1 

-74.  6 

-39.  5 

84.  0 

9344.  0 

(c)  Path  KON-ALQ 


Event:  LX+GTURP+1  34  ; 02/19/73  ; 40.  2N  , 33.  9E  ; mb=  4.  9 

Station 

Latitude 

°N 

Longitude 

°E 

* 

AZ 

4° 

4 km 

KON 

59.  7 

9.  59 

-30.  0 

24.  7 

2747. 1 

ALQ 

34.  9 

-106.  5 

-31.7 

96.  6 

10749.  7 

* AZ: 

4 

**  A : 


Source  azimuth  + for  NE 

Epicenter  distance  in  degrees  and  km 
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• KON  to  ALQ  (VLPE  station  at  Albuquerque,  New  Mexico)  ob- 

tainable from  the  event  LX+ GTURP+ 1 34. 

Again,  the  two  stations  involved  are  not  exactly  on  the  same  event- station 
great  circle  path;  however,  the  difference  is  small  as  compared  to  that  in  the 
previous  TLO-ZLP  path;  the  source  azimuths  of  the  two  stations  differing  by 
0.  4 degree  for  KON-TLO  path,  1.  6 degrees  for  KON-OGD  path,  and  1. 7 
degrees  for  KON  to  ALQ  path.  The  calculated  energy  attenuation  coefficients 
for  these  three  paths  are  shown  in  Figure  III-4. 

Among  these  four  travel  paths,  the  KON-OGD  path  lies  mostly 
in  the  Atlantic  Ocean;  the  KON-TLO  path  lies  mostly  along  the  continental 
shelf  of  the  west  coast  of  Europe;  the  KON-ALQ  and  the  TLO-ZLP  paths  are 
mixed,  more  continental  for  the  former  and  more  oceanic  for  the  latter. 
Referring  to  Figures  III-2  and  IU-4,  it  is  noticed  that  the  calculated  attenuation 
curves  are  quite  irregular  and  comparison  seems  hard.  However,  the  calcu- 
lated attenuation  coefficients  here  seem  to  be  reasonable  in  the  sense  that  they 

correspond  to  the  realistic  Q values  between  100  and  1000,  where  Q = Z7tf/k  U 

£ 

and  U is  the  group  velocity  at  frequency  f.  One  common  feature  for  these 
four  paths  is  that  the  attenuation  coefficients  for  the  10-second  period  here 
appears  to  be  too  low  relative  to  the  other  periods.  This  is  judged  from  the 
existing  knowledge  that  for  the  continental  path  the  attenuation  should  increase 
significantly  with  decreasing  period  in  the  shorter  periods  (Tryggvason,  1965); 
and  the  oceanic  path  usually  has  higher  attenuation  than  the  continental  path, 
especially  for  the  shorter  periods.  The  exact  cause  for  this  relatively  low 
attenuation  coefficient  at  10  seconds  is  not  known;  but  may  possibly  be  caused 
by  the  improper  instrument  response  correction  for  the  10 -second  period,  as 
commonly  observed  at  the  VLPE  stations  and  by  the  fact  that  the  observed  10- 
second  energy  was  relatively  small  and  the  extraction  of  the  10 -second  spec- 
trum could  not  be  made  accurately. 
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For  the  KON-OGD  path,  the  computed  attenuation  coefficient 

-5  -1 

at  40  seconds  is  unreasonably  low  (only  10  km  ) and  those  at  45  and  50 
seconds  are  negative,  physically  impossible.  One  possible  explanation  is 
that  these  two  stations  are  not  exactly  on  the  same  event-station  great  circle 
path,  the  great  circle  path  being  very  close  to  the  nodal  direction  of  the  Ray- 
leigh wave  radiation  pattern.  For  the  periods  between  15  and  35  seconds,  the 
level  of  the  attenuation  curve  obtained  here  is  compatible  with  that  obtained 
by  Tsai  and  Aki  (1969)  for  some  other  oceanic  path  across  the  Atlantic 
Ocean,  although  the  detailed  variations  are  quite  different. 

For  the  KON-TLO  path,  the  level  of  the  computed  attenuation 
curve  is  higher  than  that  obtained  for  the  other  three  paths  here.  This  seems 
unusual,  since  the  KON-TLO  path  is  mostly  continental  while  the  other  three 
are  either  entirely  oceanic  or  mixed.  This  situation,  however,  will  also  be 
observed  for  the  average  travel  path  in  the  northern  part  of  Eurasia  as  we 
apply  the  isotropic-source  method  to  the  explosion  data  in  the  following  sub- 
section. 

2.  Explosion  Data 

Long  period  Rayleigh  waves  of  presumed  underground  explosions 
from  the  U.  S.  Nevada  test  site  (NTS)  and  from  the  Russian  eastern  Kazakh 
(EKZ)  taking  place  during  the  first  six  months  of  1975  when  the  recordings  at 
three  array  stations  (ALPA,  LASA,  and  NORSAR)  were  all  readily  available 
on-line  to  SDAC  were  examined  in  order  to  estimate  the  average  Rayleigh  wave 
attenuation  for  the  travel  paths  from  NTS  and  EKZ  to  three  array  stations. 

Five  NTS  and  three  EKZ  events  were  available  and  usable  for  this  purpose. 

For  NTS  events,  the  observed  Rayleigh  wave  data  were  available  and  good  at 
all  three  array  stations;  while  for  EKZ  events,  data  quality  was  not  as  good  as 
that  of  the  NTS  and  two  events  did  not  have  LASA  recordings. 

Table  IU-3  presents  the  information  of  events  and  stations  per- 
tinent to  the  calculation.  The  involved  event-station  great  circle  paths  are 
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INFORMATION  OF  EVENTS  AND  STATIONS  FOR  NTS  AND  EKZ  EVENTS 


* AZ:  Source  azimuth  + for  NE 

**  A:  Epicenter  distance  in  degrees  and  km 


I 


plotted  in  Figures  III-5  and  III-6  for  NTS  and  EKZ  events,  respectively.  It 
is  noticed  that  these  travel  paths  are  mostly  continental  --  over  the  northern 
part  of  North  America  for  NTS  events  and  over  Eurasia  for  EKZ  events.  This 
uniformity  among  the  travel  paths  is  desirable  in  using  the  isotropic-source 
method. 


We  first  applied  ISM-A  to  both  NTS  and  EKZ  events.  The 

average  energy  attenuation  coefficients,  k , for  periods  from  10  to  50  sec- 

£ 

onds  computed  by  equations  (III-6)  and  (III-7)  are  plotted  in  Figure  III-7. 

Notice  that  the  estimated  k^  values  from  EKZ  events  are  more  scattered 
than  those  from  NTS  events.  This  is  thought  to  be  reasonable,  since  for  EKZ 
events  data  quality  was  not  as  good  as  that  of  NTS  because  the  number  of  obser- 
vation stations  were  fewer  and  the  k values  were  averaged  over  a fewer  num- 
ber  of  events.  We  next  applied  ISM-B  to  the  same  events  listed  in  Table  III- 3. 
The  estimated  k values  are  again  plotted  in  Figure  HI-7.  Referring  to 
Figure  HI-7,  it  is  fair  to  say  that  the  ISM-A  and  the  ISM-B  have  yielded  com- 
patible results.  However,  the  ISM-B  estimated  attenuation  curves  appear  to 
be  smoothed  and  to  offer  better  estimates  for  the  general  trend  of  the  attenu- 
ation curve  in  the  low  frequency  range  of  20  to  50  seconds.  For  comparison 
purposes,  the  Tryggvason  attenuation  curve  (Tryggvason,  1965)  is  also  shown 
in  Figure  IH-7. 

Again  referring  to  Figure  HI-7,  but  this  time  focusing  attention 
on  the  attenuation  curves  obtained  from  NTS  events  (NTS  curve)  and  those  ob- 
tained from  EKZ  events  (EKZ  curve),  it  is  not  hard  to  see  that  the  EKZ  curve 
and  the  NTS  curve  are  quite  different  in  two  respects.  First,  the  level  of  the 
EKZ  curve  is  much  higher  than  that  of  the  NTS  curve  in  the  low  frequency 
range  of  20  to  50  seconds.  Second,  the  EKZ  curve  has  a maximum  value  at 
40  seconds  and  decreases  toward  both  longer  and  shorter  periods  (down  to  15 
seconds);  while  the  NTS  curve  has  a minimum  value  at  35  seconds  and  in- 
creases toward  both  longer  and  shorter  periods. 
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EVENT-STATION  TRAVEL  PATHS:  NTS  TO  THREE  ARRAY  STATIONS 
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EVENT-STATION  TRAVEL  PATHS:  EKZ  TO  THREE  ARRAY  STATIONS 


(ISM-A) 


The  frequency  variation  of  the  EKZ  curve  may  not  be  reliable 
and  probably  is  not  too  significant,  since  the  estimated  k^,  values  from  EKZ 
events  are  scattered,  as  discussed  earlier,  and  the  attenuation  curve  here 
has  been  obtained  by  just  connecting  the  k values  for  nine  discrete  frequen- 
cies.  However,  the  large  difference  between  the  EKZ  curve  and  the  NTS  curve 
in  the  level  of  the  attenuation  curve  is  thought  to  be  believable,  since  it  merely 
reflects  the  large  difference  between  the  observed  signal  amplitudes  of  EKZ 
events  and  those  of  NTS  events,  as  shown  in  the  following  figure.  Figure  HI-8 
shows  the  observed  Rayleigh  wave  traces  recorded  at  NORSAR  for  4/30/75 
NTS  event  (NTS/430/75N0)  and  for  4/27/75  EKZ  event  (EKZ/427/75N0).  Also 
given  in  the  figure  are  the  values  of  the  bodywave  magnitude  (mb),  the  epicenter 
distance  (A),  and  the  maximum  recorded  amplitude  (zero-to-peak)  in  computer 
count  (CC).  It  is  noticed  that  EKZ/427/75N0  had  a half-magnitude  larger  mfe 
and  about  one-half  time  smaller  A than  NTS/430/75  (5.6  versus  5.  2 and  38 
versus  73°);  yet  two  events  produced  almost  the  same  observed  signal  ampli- 
tude (55.  2 CC  versus  51.  5 CC).  A simple  explanation  is  that  the  energy  atten- 
uation along  the  path  from  EKZ  to  NORSAR  is  much  larger  than  that  along  the 
path  from  NTS  to  NORSAR. 

Comparing  the  NTS  curve  with  the  Tryggvason  curve,  it  is 
found  that  general  trends  of  two  curves  are  quite  similar,  especially  when  a 
straight  line  is  fitted  through  the  NTS  curve  from  50  to  20  seconds.  The  gen- 
eral trend  is  that  the  k_  has  a local  maximum  of  50  seconds,  decreases 

£ 

slowly  toward  20  seconds,  and  increases  linearly  from  20  to  10  seconds. 
Nevertheless,  a small  difference  does  exist  between  the  two  curves;  namely, 
the  NTS  curve  has  a little  bit  higher  value  at  20  seconds  and  increases  faster 
toward  10  seconds. 
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CONCLUSIONS 


Rayleigh  wave  energy  attenuation  coefficients  in  the  period 
range  of  10  to  50  seconds  have  been  estimated  for  several  travel  paths  by 
applying  the  two-station  method  and  the  isotropic-source  method  to  the  earth- 
quake and  explosion  data.  In  general,  the  results  are  thought  to  be  reasonably 
realistic.  However,  no  attempt  has  been  made  to  interpret  the  obtained  atten- 
uation curves  in  terms  of  the  Q -values  of  the  earth  structure.  The  attenuation 
coefficients  computed  here  are  intended  to  be  used  for  the  travel  path  attenu- 
ation correction  in  the  far-field  source  parameter  estimation  and  possibly  in 
the  calculation  of  surface  wave  magnitude.  Based  on  the  knowledge  gained 
from  the  travel  path  attenuation  study  presented  here,  the  following  remarks 
can  be  stated: 

• The  two-station  method  is  simple;  however,  its  practical  appli- 
cation is  thought  to  be  limited  due  to  the  fact  that  it  is  hard  to 
have  two  stations  lie  exactly  on  the  same  event-station  great 
circle  path.  Serious  error  in  the  computed  energy  attenuation 
can  be  made  when  the  source  azimuths  of  two  stations  are  not 
exactly  the  same  and  are  close  to  the  nodal  direction  of  the 
surface  wave  radiation  pattern. 

• A noticeable  difference  in  the  frequency  variation  of  the  com- 
puted attenuation  curve  can  be  expected  among  different  travel 
paths.  Hence,  it  seems  that  the  attenuation  coefficients  ob- 
tained for  a specific  travel  path  by  the  two-station  method  can, 
at  best,  be  used  along  that  path. 

• The  ISM  - B applied  to  a small  number  of  observation  stations 
with  many  available  events  can  be  expected  to  yield  compatible 
results  as  the  ISM-A  applied  to  a large  number  of  observation 
stations  with  one  single  event.  When  only  the  small  number  of 
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observation  stations  with  many  events  is  available,  the  ISM-B 
can  be  expected  to  do  better.  However,  with  the  ISM-B  the 
data  quality  should  be  good  and  consistent  among  events;  other- 
wise, the  computed  k values  can  be  negative  and  physically 
unrealistic. 

The  average  Rayleigh  wave  attenuation  coefficients  estimated 
from  the  NTS  events  for  the  travel  paths  from  NTS  to  three 
arraY  stations  agree  quite  well  with  the  Tryggvason  attenuation 
curve,  although  the  former  has  a little  bit  higher  value  for  the 
short  periods. 

Based  on  the  observed  Rayleigh  wave  data  used  here,  the  large 
energy  attenuation  coefficients  obtained  for  the  path  from  EKZ 
to  three  array  stations  are  believable,  although  this  seems 
unusual  and  may  not  be  a general  situation. 
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APPENDIX  A 


Observed  surface  waves  for  the  event  L.X+SINK+S018  in  the 
1974  Sinkiang,  China, earthquake  swarm  are  shown  in  the  following  figures. 
These  figures  are  thought  to  be  representative  for  the  events  in  this  swarm. 
In  these  figures,  for  each  observation  station,  the  recorded  surface  wave 
trace  is  given  first,  followed  by  the  average  group  velocity  curve  for  that 
travel  path  and  the  uncorrected  amplitude  spectrum  which  are  obtained 
from  the  narrowband  filtering  process. 
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FIGURE  A-ld 

LOVE  WAVE  TRACE:  RECORDED  AT  TLO 
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LOVE  WAVE  TRACE:  RECORDED  AT  KIP 
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FIGURE  A-3f 

UNCORRECTED  LOVE  WAVE  SPECTRUM:  ALPA 
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RAYLEIGH  WAVE  TRACE:  RECORDED  AT  NORSAR 


FIGURE  A-4c 

UNCORRECTED  RAYLEIGH  WAVE  SPECTRUM:  NORSAR 
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LOVE  WAVE  TRACE:  RECORDED  AT  NORSAR 
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APPENDIX  B 

The  following  figures  present  the  distribution  of  the  minimum 
residuals  of  each  source  parameter  for  the  second  group  of  events  discussed 
in  Subsection  1I-C  with  their  spatial  fits. 
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